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We propose a scheme to achieve a uniform cross phase modulation (XPM) for two nonclassical light 
pulses and study its application for quantum non-demolition measurements of the photon number 
in a pulse and for controlled phase gates in quantum information. We analyze the scheme by 
quantizing a common phenomenological model for classical XPM. Our analysis first treats the ideal 
case of equal cross-phase modulation and pure unitary dynamics. This establishes the groundwork 
for more complicated studies of non-unitary dynamics and difference in phase shifts between the 
two pulses where decohering effects severely affect the performance of the scheme. 



PACS numbers: 

I. INTRODUCTION 

Optical cross phase modulation (XPM) is a specific 
variant of nonlinear optical phenomena of Kerr type in 
which the refractive index n\ of light in pulse 1 varies lin- 
early with the intensity I2 of another light field, so that 
ei = 1 + Xi + X3^2j with ei = n\ the permittivity of the 
medium, xi the linear, and X3 the nonlinear XPM sus- 
ceptibility, respectively. In ordinary optical media this 
effect is small and requires large intensities, but the pro- 
posal by Schmidt and Imamoglu pQ to generate giant 
nonlinearities using electromagnetically induced trans- 
parency (EIT) 2-4 has made very large Kerr coefficients 
possible [TJ EHS] and may even lead to nonlinear effects 
at the single-photon level [5HT5]. 

XPM is a strong candidate for the design of opti- 
cal quantum controlled-phase gates (CPG) for photonic 
quantum information processing |16H18| . in which two 
single-photon pulses would become entangled. For suffi- 
ciently high fidelity, XPM would allow the construction of 
deterministic gates, as opposed to the non-deterministic 
optical CPG that are based on linear optics [H] and es- 
tablish entanglement through measurements [2U]. It has 
been suggested that double EIT - EIT for both pulses 
with matched group velocities - would generate the max- 
imal XPM phase shift because matched group velocities 
for both photons would maximize their interaction time 
[SHU]- However, XPM based on double EIT still faces 
some challenges: (i) to achieve sufficiently high intensi- 
ties at the two-photon level, the photon pulses must be 
tightly confined in the transversal direction; (ii) if the 
nonlinear medium has a finite response time, the matter- 
light interaction unavoidably induces noise. It was shown 
by Shapiro [21 that for a large class of XPM models this 
would limit the fidelity of a CPG to only about 65%. 
Further theoretical studies on this topic confirmed these 
results [2"2Tl25| . (iii) For matched co-propagating pulses, 
each point of the pulse will experience a different XPM 
phase shift because the intensity, and hence the refractive 
index n\, varies over the shape of the pulse. This would 
severely affect the entanglement between the pulses. 




FIG. 1: Scheme to achieve uniform XPM: two light pulses 
, tp2 travel at different group velocities v± > V2 through an 
XPM medium such that ^\ can overtake 1^2- d denotes the 
initial distance between the two pulses. 

In this paper we will address the third problem and do 
not deal with problem (i) and (ii), i.e., we will assume 
that transverse confinement of the photon pulses can be 
achieved by some means, such as hollow core fibers |26l - 
|2T)] or nano wires [3D], and that the medium's response 
time is so short that the polarizability of the medium 
reacts instantaneously to a photon's electric field. The 
omission of the noise associated with a finite response 
time is also made for clarity because, despite that we are 
working in the instantaneous regime, we will show that 
very similar effects may appear if the two light pulses 
have different group velocities. To solve problem (iii) 
we extend an idea of Rothenberg |31j who demonstrated 
a uniform XPM phase shift for classical pulses in bire- 
fringent fibers: two light pulses travel through a XPM 
medium with different group velocities v\ > i>2, see Fig. [I] 
Pulse 2 (blue) reaches the Kerr medium first but pulse 1 
(green) is faster and leaves it first. While the pulses over- 
lap they interact via the Kerr effect, which is proportional 
to the pulse intensity. Because pulse 1 overtakes pulse 2 
the acquired phase shift will be averaged over the pulse 
shape and thus be nearly uniform. 

Here we generalize this idea to characterize the propa- 
gation of quantized light pulses and show that a uniform 
phase shift can also be achieved for quantum interference 
effects. We will show that quantizing the classical XPM 
equations is generally a subtle problem which may make 
the introduction of noise terms similar to those in prob- 
lem (i) necessary. We exploit our results to suggest exper- 
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iments for quantum non-demolition (QND) measurement 
of photon numbers and CPG. 

This paper is organized as follows. In Sec. [n] we will 
present a phenomenological quantum model for XPM be- 
tween two pulses with different group velocities and dis- 
cuss the two cases of unitary and non-unitary dynamics 



in this model. In Sec. Ill we present proposals for QND 
measurement of photon numbers and CPG and show that 
they can have a high fidelity in this case. The analysis of 
these proposals for the non-unitary case is presented in 



Sec. IV Several appendices contain details of our deriva- 
tion. 



II. PULSE PROPAGATION WITH 
MISMATCHED GROUP VELOCITIES 

We consider two quantized light pulses that travel 
through an XPM medium with different group velocities 
vi > v 2 - We assume that the light pulses are transver- 
sally confined so that we can restrict the model to one 
spatial dimension. In classical fiber optics this situation 
is commonly described by the set of equations [32TI34] 
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with imaginary unit i = \J — 1. Here £i is the field am- 
plitude of pulse i = 1,2, = \£i\ 2 is proportional to the 
intensity, and ji is the self phase modulation (SPM) coef- 
ficient. The XPM coefficient is given by 2ji. This model 
is local and instantaneous; i.e., the light field £i(z,t) is 
only affected by £2(2, t) and vice versa. 

In the context of SPM, Joneckis and Shapiro have 
shown that in a local and instantaneous quantum model, 
where £i represents a field operator, infinite zero-point 
fluctuations do occur [35]. They suggested taking a fi- 
nite response time of the medium into account to avoid 
singularities. We instead propose a model in which the 
medium's response is still instantaneous but spatially 
nonlocal, i.e., £i(z,t) can be affected by £2(2', t). The 
most general macroscopic model that describes such an 
instantaneous interaction between two light pulses of dif- 
ferent group velocities is then given by 
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d t + d z ) £i(z,t) 



dz' h-i{z' ,t) 



x Vi{z-z!)£i(z,t). 



(2) 



Here and henceforth £i denotes the field operator and 
Ii(z,t) = £j(z' ,t)£i(z' ,t) the intensity operator. The 
quantities Vi(z-z') are spatially nonlocal interaction po- 
tentials between the two pulses generated by the atomic 
medium. To avoid causality violation its support must 
be much smaller than the wavelength of light. We will 
assume that Vi{z — z') is essentially zero if the distance 
between z and z' is much larger than the Bohr radius. 



One can consider the introduction of nonlocal potentials 
as a regularization of the classical theory ([I]) by smearing 
out the point interaction. In the limit of sharp potentials, 
Vi(z) — ViS(z), the regularized model reduces to Eq. ([I]) 
for XPM coefficient 7, = Vi/2. 

Multi-photon pulse dynamics should account for SPM, 
but we omit this effect here in order to focus on XPM. 
When XPM is placed in the context of interferometry, 
which gives phase an operational meaning, pairs of non- 
linear media in both paths can offset SPM |36| . Alter- 
natively SPM is not present for single-photon pulses be- 
cause a single photon cannot induce SPM on itself. 

The multi-mode field operators £i(z) satisfy the com- 
mutation relations 
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5ijS(z - z') = r-jSijS^z - z'), 



(3) 



with uji the central frequency of pulse i and A the trans- 
verse area of the pulses. The field is assumed to be in 
a single mode with the i th longitudinal mode function 
given by tpi (z) . Photons in this mode are annihilated by 
the operator 
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Vv J- 



dz ip* (z — u) £i{z, 0) 



(4) 



For later use we have included a possible shift u of the 
wavefunction. If the shift is zero we will sometimes sup- 
press this notation and write = [0] . 

If Eq. ([2| is taken to be the dynamical equation of the 
quantum fields £i one has to be careful with its interpre- 
tation. It can be shown that Eq. ^ can be expressed as 
dt£i = — i[£i,H] for some Hamiltonian H if and only if 
v\V\{z) = v 2 V2(z). In other words, a unitary evolution 
of the quantum fields will only occur if the two interac- 
tion potentials are proportional to each other, with the 
proportionality factor given by the ratio of the group ve- 
locities. In absence of absorption or other decohering 
processes one therefore would assume that in the quan- 
tum model this relation must be fulfilled. However, in 
the corresponding classical models for XPM pulse prop- 
agation this is generally not the case [32H34j . 

In our subsequent analysis we will therefore first as- 
sume that the evolution is unitary and analyze our pro- 
posals within this framework. In Sec. |IV| we then will 
develop a consistent phenomenological quantum model 
for non-unitary evolution using techniques of open quan- 
tum systems and reassess the proposals for this case |37| . 



III. XPM FOR UNITARY EVOLUTION 

To devise a scheme for CPG and QND measurements 
of the photon number we will assume that the dynamics 
is unitary; in this case the interaction potentials can gen- 
erally be written as Vi{z) = V(z)/vi, with V(z) a Hamil- 
tonian interaction potential. However, we will make use 
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of this relation only after deriving general results in order 
to facilitate the extension to non-unitary dynamics. 

The solution of Eq. ([2| for the propagation of light in 
an infinitely extended XPM medium is 



£i(z,t) = U 3 -i(z,t)£i(z - Vit,0), 
Ui(z,t) =exp ( -i / h(z',0) 



(5) 



•«'3- 



Vi - V 3 -i 



x 6 3 -i (z- z 1 - v 3 -it , [vi - v 3 -i]t) dz' \ , 

(6) 

%(z 1 ,z 2 )= ( 2 dz' Vi{ Zl - z') . (7) 



The function 9j,(zi,z 2 ) will play a central role in deter- 
mining the Kerr phase shift. One important property is 
that for a symmetric potential, Vi(—z) = Vi(z), we have 
9i(-zi, -z 2 ) = -d i (z 1 ,z 2 ). 



A. Quantum non-demolition measurement of 
photon numbers 

As a first application we consider a QND measurement 
of the number of photons [35] in mode 2. This can be 
accomplished by sending a strong classical pulse in mode 
1, which can be described by a coherent state \a), to- 
gether with an n photon pulse \n) in mode 2 through the 
XPM medium. Single-mode treatments predict that in 
this case the phase of the classical pulse will be shifted 
by an amount that is proportional to n. Here we show 
that the uniform cross phase shift that we suggest will 
accomplish precisely this. 

The initial state of the two light pulses before they 
start to interact takes the form — \a) ® |n) where 
the two kets refer to mode 1 and 2, repectively. Using 
the shift operator D\(a) = exp(aa{ — a*a\) this can be 
expressed as 



(8) 



The complex amplitude of the classical pulse at time t 
and position z is given by (V>|£i(z, t)\ip). Its phase can 
be measured using homodyne detection |39j . for instance, 
which more specifically measures th e ob servable X($) = 



e m Si(z, t) + e- M £l(z, t). Using Eq. (IA3J it is easy to see 



that 



D\ {a)£i 0, 0)L>i (a) = £ i (z, 0) + a^r/Vi (*) 
Exploiting this and solution ^ we get 



(9) 



n 



, ^i(*-«i*)<0|a5(55[0;«,i]) n |0), 



(10) 



where we have used Eq. (A4| and introduced an XPM- 



modified annihilation operator 
1 [°° 

ai[u; z , t] = — — / dz ip* (z — u) £ x [z , 0) 

V V J—oo 



x e 



-i(z— z' -V 3 _it,(Vi— V 3 _i)t) 



(11) 



The physical interpretation of Eq. (Ill is that the 



wavepacket ip*(z' — u) is multiplied by a spatially varying 
phase factor that is given by the exponential in Eq. ( 11 ) 



and incorporates the effect of XPM on the light pulses. 



The expectation value for field 2 in Eq. ( 10 ) can then 
be reduced to 



a^(alp,z,t}) n \0)=n\[a 2 ,al[0-,z,t]} n , (12) 



Wan nil 



so that 



{£\ (z, t)) = a^fqi^i (z - vit) 

x (^J dz' \Mz')\ 2 e iv ^ ei(z '~ z+Vlt ' Avt) J\ 



(13) 



with Av = v\ — v 2 . 

To better understand the physical implications of 
Eq. ( [i~3| w e consider the specific configuration depicted 
in Fig. [1] The classical mode tpx is initially centered 
around the origin while the n-photon pulse in mode ?p 2 is 
a distance d to the right of it [H] . Both pulses are mov- 
ing to the right, but pulse 1 is faster. The first line in 
Eq. ( 13 ) is basically the amplitude of the classical pulse 
ipi at time t in absence of the XPM medium. The pulse 
is centered around z — v\t. Hence, to achieve an maxi- 
mum phase contrast, we should observe the field at this 
point. The exponential in Eq. (13) is then proportional 
to O^z^Avt). 

In Appendix |B| it is shown that for any potential Vi(z) 
that is consistent with causality, the function 9\{z' , Avt) 
is nearly constant between the lines Av t = z' and z' = 
and zero outside of this range. The support of the ini- 
tial wavepacket ip 2 (z') is in the area z' > and peaked 
around z' = d. For sufficiently large times, such that 
Avt ^> d, the support of ip 2 (z') is therefore completely 
inside the domain where 6\{z' ,Avt) ~ Vi, with the con- 
stant 



dz Vi(z). 



(14) 



This condition corresponds to the requirement that the 
classical pulse 1 had enough time to overtake the n- 
photon pulse 2. Because the mode function ip 2 is nor- 
malized we thus find 



(£i{vit,t)) 



:a0?Vi(O)e- in ^ 
Av 



(15) 
(16) 
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FIG. 2: Construction of a controlled phase gate for two time- 
bin qubits using uniform XPM. 



We remark that (f>\ — 02 for unitary dynamics. Eq. (15) 



is precisely the result that one would obtain in a single- 
mode treatment. Hence, for unitary dynamics a mis- 
match in the group velocities of two pulses that allows 
one pulse to overtake the other in a Kerr medium will 
result in a uniform phase shift for QND measurements of 
the photon number. In practice, overtaking will impose a 
minimum requirement on the length of the Kerr medium 
and SPM will lead to a distortion of the signal. 



for (6zP,5z!j 2) ) = {0,0),(0,Sz),(Sz,0),(dz,Sz) for 
(ij) = (0,0), (0,1), (1,0), (1,1), respectively. The fac- 
tors vit represent an explicit time dependence in the def- 
inition of the logical states such that they are co-moving 
with the wavepackets at group velocity u,. 

We consider the case that the two photons are ini- 
tially prepared in a pure two-qubit state \ip(0)) — 
J2kl c ki\kl(0))- To characterize the performance of the 
CPG we need to obtain the logical components Pij-u (t) — 
(ij(t)\p(t)\kl(t)) of the density matrix p(t) at time t. 
A somewhat tedious calculation that can be found in 
App. |C 1| leads to 



with 



X, 



x e 



Pij-ki{i) = X* kl Xij, (19) 



dzi dz 2 |t/>i(zi) ip2{z2 + d)f 



-ir)^e 2 (z2~z 1 +d+Szf ) -Szl 11 ,Avt) ^Q) 



B. Controlled-Phase Gate 

Our proposal to build a CPG extends previous designs 
to build quantum gates using XPM 16, 40] by ensuring a 
uniform XPM phase shift. In a controlled phase gate, two 
qubits with logical basis states |0), |1) are manipulated in 
such a way that the state acquires a phase shift j3 only 
if both qubits are in the logical state |1) ® |0). In other 
words, a CPG is a unitary operator that maps the basis 
states \kl) = \k) <x> \l) to [12"] 



The phase factor in Eq. (20) determines the XPM ef- 



U C p G \kl)=e ik ^- l ^\kl) , k,l = 0,1 



(17) 



With the uniform cross phase shifter discussed above, 
such a gate can be implemented for time-bin encoded 
photonic qubits as indicated in Fig. [2j In this case each 
optical field £ n carries one photon. For each photon the 
logical basis states are encoded in two wavepackets, one of 
which is delayed compared to the other. The logical state 
|0) for qubit n (n — 1, 2) can be encoded in wavepackets 
ip n {z) so that photons in this state are annihilated by 
operator a n = a n [0]. The logical state |1) can be encoded 
in wavepackets that are shifted to the right, ip n {z — 5z); 
photons in this state are annihilated by the operators 
Orjtfcj. A CPG can then be constructed by arranging 
the initial distances between the wavepackets in such a 
way that the only states that do overlap inside the XPM 
medium are state |1) of the first qubit and state |0) of the 
second qubit. For all other logical states the wavepackets 
in the two field modes never overlap so that there is no 
XPM effect. 

The logical basis states for two time-bin encoded qubits 
at time t are given by the set 



feet. We now prove that it is only nontrivial for the 
logical state 1 10) because only in this case the associated 
wavepackets overtake each other. 

The phase factor is proportional to the function 62 {zi — 

z\ + d + Szj — i Avt) which, according to Appendix 
|B"} is nearly constant for a given range of its two vari- 
ables. Ignoring the variables Zj for the moment we con- 



sider 6*2 (d + Sz 



(2) 



Avt) instead. For the four 



states 00, 01, 10, 11 the first argument takes the values 
d.d + 5z, d — Sz, and d respectively. If we choose the 
time t for which the photons are interacting such that 
d - Sz < Avt <d,d + Sz then 2 = for ij = 00, 01, 11 
and 82 = V2 for ij — 10. Hence only the latter combina- 
tion will experience an XPM phase shift. This conclusion 
also holds if we re- introduce the variables z, because they 
only vary over the support of the two wavepackets, which 
is smaller than d, Sz; the value of 8 2 therefore does not 
change. Thus we arrive at 



X{j CijC ' 



and the density matrix elements become 



Pij;kl{t) 



CijC kl e 



* p -i(i(l-i)-fc(l-O)02 



(21) 



(22) 



Let us now work out the performance of the CPG 
by calculating the concurrence |43j of the density ma- 
trix (22) for a specific initial product state | -0(0)) 



(|0) +]?)) <S> (|0) + |l))/2, which corresponds to cy = \ 
for all The concurrence C of the density matrix 



\ij(t)) = a\[Szl 1} + Vl t] 4M 2) + v 2 t] |0), 



(18) 



(22) then evaluates to C = \ sin(</>2/2)|. This indicates a 
perfect CPG because the initial product state is trans- 
formed into a maximally entangled state provided one 
can achieve a large phase shift <p2 = tt. 
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IV. XPM WITH NON-HAMILTONIAN 
COUPLING 

In the previous section we have devised a scheme to 
create uniform XPM for quantized light pulses for the 
special case of unitary dynamics for which viV\(z) — 
v%V2(z). The theory behind this scheme corresponds to 
a regularized quantization of the classical theory ([I]) for 
the case that the XPM coefficients are related by 1*171 = 
i>272- However, in classical optics this relation between 
the two XPM coefficients is not generally adopted. In 
this section we therefore study the quantization of the 
classical model with U171 7^ V272 and its implications 
for quantum information processing. Because the model 
needs to be regularized, we consider the general case with 
two non-local interaction potentials that fulfill v\ V\ (z) 7^ 

V2V 2 {z). 

The discussion in Sec.[TI]lias revealed that the quantum 
dynamics in this case must be non-unitary because there 
is no Hamiltonian that can generate the equations of mo- 
tion. As Eq. |2]) is based on a phenomenological model 
that does not explicitly include absorption or dissipation, 
the origin for this non-unitarity cannot be identified in 
an unambiguous way. There might be implicit absorp- 
tion processes hidden in the model, although we strongly 
doubt that this is the case because the dynamics does 
not have a structure that is comparable to the Lindblad 
form [33]. Another possibility is that the external fields 
that usually are needed for EIT-based XPM media may 
induce energy fluctuations that lead to non-Hamiltonian 
dynamics. Yet another possibility is that the direct quan- 
tization of the classical theory ignores the averaging 
processes that are associated with a macroscopic decrip- 
tion of electrodynamics [35] . Such averaging procedures 
lead to loss of information, which may induce decoher- 
ence in a quantum system. 

Despite this ambiguity with respect to the cause of 
non-unitary dynamics it is worthwhile studying this sit- 
uation because one can gain a better understanding of the 
quantization of a new class of classical models. Further- 
more, even if the theory that we develop below cannot 
include the microscopic details behind the XPM interac- 
tion it nevertheless should help to estimate the effects of 
a mismatch between the XPM coefficients. 

In the case ViVi(z) V2V~2(z), Eq. ^ still provides 
an exact solution to the dynamical equation of motion. 
However, this solution is inconsistent with basic require- 
ments of quantum field theory. To illustrate this point 
we consider the equal-time commutation relations be- 
tween the Heisenberg field operators, which should agree 
with the commutation relations between the respective 
Schrodinger operators. In our case this principle is vio- 
lated because 



despite [£ 1 (z,0),£ 2 (z',0)] = 0. 

The conventional way to deal with non-Hamiltonian 
dynamics in the Heisenberg picture is to introduce 
Langevin noise operators. In the context of nonlinear 
optics this has first been done by Boivin et at |46j in the 
case of time-dependent self-phase modulation. 

We now show that a similar approach can also be made 
for instantaneous XPM between two light pulses with dif- 
ferent group velocities. Following Ref. |46j we introduce 
Hermitian decoherence operators rrii(z,t) such that the 
dynamical equations are modified to 

(jdt + d)j £i{z,t) = -i( j dz'I 3 -i(z',t) (25) 

xVi(z- z) + mi(z,t)^j£i(z,t) . 

The decoherence operators commute with all field op- 
erators. We adopt the frequently used assumption of 
(^-correlated decoherence, 

[ mi (z, t) , mj (z', f )] = if) - V ^ 1 ) S(t ~ t>). 

_ (26) 
The solution of Eq. ( 25 1 can be written as 

£i(z, t) = Ri(z, t) U s -i{z, t) £i{z - Vit, 0), (27) 



Ri(z, t) = exp ( — i«j J dt' m,i(z — Vi(t — t'),t') 



(28) 



Using the commutation relations ( 26 1 it is not hard to 
show that 

R\(z, t) R 2 (z / , t) R x (z, t) = R 2 {z', t) e^C*' -*•*), 
Rl(z,t)R 1 (z',t)R 2 (z,t) = R^z 1 ,t) e-^ {z - z ' ' t] . (29) 



As a consequence, solution (27 1 fulfills the equal-time 
commutation relations. Introducing the decoherence 



term in Eq. (251 has thus led to a consistent quantum 



field theory of XPM for light pulses with different group 
velocities. 



A. Controlled-Phase Gate 

To characterize the performance of a CPG in the case 
of non-unitary dynamics we can repeat the steps dis- 
cussed in Sec. Ill B with solution ^ replaced by Eq. ( 27 1 . 
Finding the logical density matrix elements must then be 
done within an open system approach that is described 



in App. C 2 As a result, the factors X^j become opera- 
tors on the Hilbert space of the environment (on which 
the decoherence operators act) and Eq. (21) is replaced 
by 



£ 2 (z',t)£ 1 (z,t)=£ 1 {z,t)£2(z',t)e- i *( z '- z > t \ (23) Xij = f dz x dz 2 |Vi(zi) M*2 + d)\ 

$(z' -z,t)=-^- {vi6 x (z' - z + Avt, Avt) 
Av 

- V202(z' - z + Avt, Avt)), (24) 



x i?i(zi + vit + 8z ( l ] ,t)R 2 {z 2 + d + v 2 t + 5zf ] ,t). 

(30) 
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The logical density matrix elements take the form 

p ij;k i(t) =Tr E (pv(0)Xl l X ij ), (31) 

with Tig indicating the trace over the environment de- 
grees of freedom and pe(0) the initial state of the envi- 
ronment. 

To keep the discussion concise we assume that the 
photon pulses are very sharp. We then can replace the 
squares of the wavefunctions by S distributions and get 

Xij = aje-^-^R^t + Szjp, t)R 2 (d + v 2 t + 5zf\t). 

(32) 

The density matrix elements then become 

PirMit) = Cij^e-^-^-^-^Nii^t), (33) 

with the Hermitian decoherence matrix 

N irM (t) = (R\{d + v 2 t + 5z\ 2 \t)R\{v l t + 8z^\t) 

x R x {v x t + , t)R 2 (d + v 2 t + 5zf ] ,t)Y 

(34) 



In App. C 3 ) we show that it can be written as 



N = 
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C2 




N 00 ., n \ 




r* 
c 2 
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Aoi;10 


Cl 




c *g-i(4>l-<fe) 


N* 

JV 01;10 
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c 2 
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N* 


r* 
c l 


r* 
'2 


1 / 



Rl(v n t + 6z,t)R 1 (v n t,t)\. 



(35) 
(36) 



Hence the influence of decoherence in our model can be 
described by four complex parameters Nqo-u, Nqi-,10, and 
c n . This is a consequence of the approximation of very 
short pulses. For extended pulses additional decoherence 
contributions are anticipated. 

We again calculate the concurrence of density matrix 
|33| for the initial product state = | for all i,j. 
The specific values of the parameters iVoo ; ii, AToi ; io, and 
c n depend on the detailed decoherence model that one 
employs so that the predictions may vary substantially 
depending on the assumptions behind the model. We 
present one particular decoherence model in Appendix 
|D| However, all decoherence models must be consistent 
with all eigenvalues of the density matrix ( 22 ) taking val- 
ues between and 1. 

To simplify the discussion we consider the special case 
4> 2 = 7r for the oft-used assumption that the two inter- 
action potentials Vi(z) are equal and independent of the 
group velocities [33J [M] so that V\ = V 2 . The phase fac- 
tor that appears in the decoherence matrix (351 is then 
4>\ — 4> 2 = irAv/v 2 > 0. For a given value of <f>i — <f> 2 
we can derive a "minimal decoherence model" by find- 
ing those values of the parameters Nqo-xx, iVoi ; io, and c n 
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FIG. 3: Concurrence as a function of the XPM phase shift 
for minimal decoherence (dots) and the decoherence model 
presented in Appendix |I5| (solid line). 



which maximize the concurrence C for a consistent den- 
sity matrix. Fig. [3] shows the concurrence for this min- 
imal decoherence model as a function of 4>i — 4>2- The 
dots are found through a random search for the param- 
eters iVoo ; ii, ^Voi;iOi an d c n of the minimal decoherence 
model. Each dot is based on a sample of typically 10 7 
random events. The error in each value is estimated to 
be about 5%. The concurrence decreases with <pi — cf> 2 
and becomes zero for cf>i — <f> 2 = w. In this case deco- 
herence completely destroys the entangling capacity of 
the XPM interaction. The numerically determined den- 
sity matrix then corresponds to a nearly equal mixture of 
two highly entangled states. Realistic decoherence mod- 
els would typically predict a less than optimal perfor- 
mance. For instance, the blue line in Fig. [3] shows C for 
the decoherence model of App. [D] under the fairly opti- 
mistic assumption Avt = 4y/irw, where w is the width of 
a Gaussian potential V\ (z) — V 2 (z) = F(0) exp(— z 2 /w 2 ). 
For more realistic values Avt 3> w |47) the concurrence 
would be non-zero only in a very narrow range around 

01 = 02- 

Our phenomenological model therefore suggests that a 
CPG may not be achievable with XPM unless <\>\ is very 
close to 02- In most proposals for XPM based on double 
EIT this could be achieved by a suitable preparation of 
the atomic gas, although it may require fine tuning of 
parameters like magnetic fields or pump field intensities. 
A better way to achieve <j)\ = <j) 2 would be to find a 
system in which this is guaranteed through microscopic 
symmetries. 



B. Quantum non-demolition measurement of 
photon numbers 

The scheme for a QND measurement of the photon 
number discussed in Sec. |III A| can easily be extended 
to the non-unitary case by calculating the measurement 
signal (£i(z,t)) using solution (27) instead of The 
only change is that the measurement result ( |15[ ) is multi- 
plied by (Ri(z,t)). This factor depends strongly on the 
specific decoherence model but generally will lead to a 
decrease in the contrast of the phase measurement. 

To give a rough estimate we consider the decoherence 
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model of App. [D] in which the noise is generated by an 
environment consisting of a harmonic oscillator field with 
a flat frequency distribution. We then have 



(Ri(z,t)) = exp 



V Vi 



4^3- 



-AvtV{0) 



(37) 



The decoherence amplitude (Ri(z,t)) can vary between 
the optimal value 1, if the interaction potential V(z) = 
Vi(z) — Vi(z) vanishes at the origin, and very small val- 
ues for a sharply peaked potential. We remark that this 
result also demonstrates the necessity of a non-local po- 
tential to avoid the diverging quantum fluctuations in the 
case V(z) oc S(z) [35] . 

It is instructive to evaluate this for a Gaussian poten- 
tial V(z) = V(0)exp(-z 2 /w 2 ) where V(0) > and the 
width w should be not substantially larger than an atom. 
We then find 4>\ — <f>2 ~ \/TrV(0)w/ri so that 



(Ri(z,t)) w exp 



Avt 



(38) 



The first factor is proportional to the XPM phase shift 
difference. The last factor Av t/w is typically much larger 
than one because the length Av t must be larger than 
the width of the light pulses (t is the time needed for 
one pulse to overtake the other). Hence, for very short 
(femtosecond) pulses the exponential may be of the same 
size as the XPM phase shift. Typically, however, it will 
be much larger. We therefore conjecture that for the 
QND measurement to be successful the XPM phase shifts 
<fii in the two light pulses must also be nearly equal. 



V. CONCLUSION 

In this paper we have devised a scheme to generate a 
uniform XPM for quantum information processing and 
applied it to characterize the fidelity of a QND mea- 
surement of the number of photons in a light pulse and 
the performance of a deterministic CPG. The analysis 
is based on a regularized quantization |2]) of a common 
phenomenological model for the classical XPM effect be- 
tween two light pulses with different group velocities. 

When the pulses overtake each other while traveling 
through a nonlinear medium the induced XPM phase 
shift (pi, i = 1, 2, is uniform over each pulse. For (fix = <p2 
the dynamics is unitary and predicts perfect fidelities for 
QND measurement and CPG. For <fii ^ (j>2 a decoherence 
effect must be introduced to keep the model consistent, 
and this decoherence significantly affects the fidelity of 
the gate and the QND measurement. 
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Appendix A: Useful commutation relations 

In this appendix we will summarize a number of com- 
mutation relations that will be useful in deriving our 
main results. 

The effect of XPM on the pulse propagation is gov- 
erned by the action of the unitary operators £/j = Ui{z, t) 
of Eq. ^ on the field operators, 



U}£i(z',0)Ui =f i (^ / ,0)exp( 



irj- 



V3- 



x 3 -i(z - z' - v 3 -it, (vi - w 3 -i)*)J , 



(Al) 



U&izf, 0)UJ =8i{z',0) exp (ir) 



V3- 



Vi - v a -i 



x 6 3 -i(z - z' - v 3 -it, (vi - v 3 ^i)t)j . 



(A2) 

The effect of XPM on the annihilation operator for 
photons in mode i can be derived from the commutation 
relations 



[<H M ,£J(z,0)] = yfr\i)* (z-u), 



(A3) 



and 



ai[u]Ui(z,t) = a,i[u;z,t], (A4) 
with <ii[u]z,t] of Eq. (111]). 



Appendix B: Simplification of the phase shift 

If the potential Vj,(z) has a finite range r, so that it 
vanishes for \z\ > r, then 9%{z\,Z2) has generally a shape 
similar to Fig. [4] It is essentially zero outside the "edges" 
depicted in the figure and has an extended plateau be- 
tween the lines Z\ — and Z\ — z 2 . The value on the 



plateau is given by sign(z2) Vi, with Vi defined in Eq. ( 14 1 



The only region where #1(21,22) varies significantly are 
bands of width 2r around the lines z\ = and z\ = z 2 , 
and only this transition region depends on the details of 
the potential. If the range r is much smaller than the 
optical wavelength we can ignore the transition region 
so that 9i{z\,z 2 ) can be replaced by the plateau value 
between the lines, 



9 l {z 1 ,z 2 ) » (000 9{z 2 - z x ) - 0{- Zl ) 9( Zl - z 2 )) Vi 
with 9{z) the step function. 



(Bl) 



Appendix C: Derivation of the action of a CPG on 
the density matrix 



1. Unitary dynamics 

The evolution of the density matrix p can be expressed 
with the aid of a unitary operator U(t) such that p{t) = 
U(t)p{<S)U^(t), with 

p(o)= E cu^ii'fmk'i'm. (ci) 

i' ,j\k,'l' 

The logical matrix elements of 
matrix can then be expressed as 
(ij(t)\U(t)p(0)W{t)\kl(t)). Using the definition ^ of 



the density 

Pij;kl(t)_ 



the logical basis states we are then led to consider the 
expression 

(ij(t)\U(t) = (OlC/^a!^ 1 ' + vit] a 2 [5zf ] + v 2 t] U{t). 

(C2) 

We here inserted an additional factor W(t) which does 
not change the result because U(t)\0) = |0). The op- 
erators ai[z] are defined using the Schrodinger operator 
£i(z, 0), see Eq. Q. The effect of ■ ■ ■ U(t) is now to re- 
place the Schrodinger operator by the respective Heisen- 
berg operator £i(z,t) so that 



We now concentrate on the factor 



= (0| ai [<5z 2 (1) +v 1 t,t] a 2 [5zf ] + v 2 t,t]\i>{0)) 
— J dz\ dz 2 il>*(zi — v\t — Szf^) 

x r 2 {z 2 - v 2 t - 5zf ] ) (0|C/ 2 Oi,t)^iOi - Vxt,Q) 
x «7iO 2 ,t)4(22-^i,0)|V(0)) 



V 



dzi dz 2 tp*(zi — v\t — ) 



x i>*{z 2 - v 2 t - Szf) (0\U 1 1 (z 2 ,t)£ 1 (z 1 - vit,0) 

x U^t)^ - v 2 t,0)\^(0)) 
1 



rh 



V 



x e 



dz[ dz' 2 ^t(z[)^(z' 2 + d) 

-iri^e2{z' 2 -z' 1 +d+8zf>-Sz!> 1) ,Avt) 



x <0|£f a (z[ + Szl 1 ' , 0)8 2 (4 + d + Sz { /> , 0)|V'(0)). 

(C6) 

The expectation value in the last line can be considerably 
simplified. 



(2) 



(0|fi (4 + 5z[ 1] , 0)£ 2 (4 + d + Szf , 0) |</;(0)) 

= E c fc ;<0|fi(4 + 0)4(4 + d + Szf,0) 



(2) 



k,l 

xa\[Sz^]al[5zP]\0) 

k,l 

x ^ 2 (4 + d + ,5zf ) -8z\ 2) ). 



(C7) 



The wavefunctions in Eq. (C6) imply that z[ <C 6z be- 



cause 4 can only vary over the support of the wavefunc- 



tions. This in turn implies in Eq. (C7) that Sz 



and Szj 2 ^ = Sz^' because otherwise the argument of the 
wavefunctions would be outside of their support. Hence, 



5z 



(i) 



(0|£i (4 + 5zY> , 0)£ 2 (z' 2 +d + 5zf ] , 0)| V(0)) 

= 77^1 (4)^2(4 + <*)> (C8) 



from which Eq. ( 20 ) follows 



(ij(t)\U{t) = (0| ai [^ (1) +v 1 t,t]a 2 [Szf ) +v 2 t,t], (C3) 
1 [°° 

ai[z,t] = — — I dzij}*(z — u)£i(z,t). (C4) 

V*j J-oo 

We thus find for the logical matrix elements 

Pir,kl{t) = (^(0)\al[Sz ( l 2) +v 2 t,t]a\[Sz ( k 1) + Vl t,t] |0) 

x (0\a 1 [Sz ( i 1) + Vl t,t] a 2 [5zf ] + v 2 t,t]\ip(0)). 

(C5) 



2. Non-Hamiltonian Dynamics 

The derivation of the logical density matrix elements 
for non-unitary evolution follows essentially the same 
steps as in the unitary case described in App. |C 1| How- 
ever, because the decoherence operators act on a differ- 
ent Hilbert space He associated with the environment, 
the calculations have to be made on an enlarged Hilbert 
space Hs+e — Hs ® %e that includes our system of pho- 
ton pulses as well as the environment. The dynamics 
on this enlarged Hilbert space is described by a unitary 
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operator U(t) and we assume that initially system and 
environment are unco rrected, ps+e(0) = p(0) ® Pe(0), 
with p(Q) of Eq. (CI I. The logical matrix elements are 
then given by 

Pir,ki(t) = T* E ({ij(t)\U(t)p(0)®p E (0)tf(t)\kl(t))), 

(C9) 

with Ti'e the trace over the environment. 



We again start the derivation with term ( C2 ) , but now 



U refers to the unitary evolution on the total Hilbert 
space. If U leaves the ground state invariant, in the sense 
that U\0) ® \E) = |0) ® U E \E) for any state \E) of the 
environment and a fixed unitary map U E that acts on the 
environment, then we can again make the transition to 
Heisenberg operators as in Eq. ( C3 ) . The logical matrix 
elements then can be written in the form (31) with the 
environment operators 

Xij = (0|ai[^ (1) +v 1 t,t] a 2 [6zf +v 2 t,t]\i>(0)) 

= ^ J dz\ dz 2 Ri(zi,t)R 2 (z2,t)ipl(z 1 - v x t - 5zf^) 
J 2 ) 



x r 2 (z2 - v 2 t - 5zf ] ) e -iv^e 2 (z 2 -z 1+ Avt,Avt) 

X (0|£i(3i - Vxt,0)8 2 (z 2 -V 2 t,0)\^(Q)). 



(CIO) 



This expression corresponds to Eq. ( C6 ) . The argument 



that leads to Eq. ( C8 ) can be repeated, resulting in 



X; 



dzi dz 2 \ipx(zx) ip 2 (z 2 + d)\" 



x R^zi +vxt + 5z\ L \t)R 2 {z 2 + d + v 2 t + 5z) z >,t) 



x e 



-ii 1 ^e2(z2-z 1 +d+8zf ) -8z[ 1 \Avt) 



(2) 

(Cll) 



which corresponds to Eq. (20 1. The discussion of the 



phase factor that leads to Eq. (21 1 is unaffected by the 
decoherence terms so that Eq. ( 30 i can be deduced. 



3. Derivation of the decoherence matrix 

The decoherence matrix can be evaluated without 
further assumptions in a number of cases. Because 
R\ (z , t) Ri(z , t) = 1 it is easy to see that the diagonal 
elements are equal to one. This also guarantees trace 
preservation. For the same reason, whenever i = k the 
decoherence operators simplify to 



N m i{t) = (^Rl(d + v 2 t + 8z\ 2 \t)R 2 {d + v 2 t + 5zf\t) 

(C12) 



Furthermore, for I — j R\ and R 2 have the same argu- 
ment so that we can make use of Eq. ( 29 1 : 



N ij;kj (t) = {R^id + v 2 t + 8zf\t)R\{ Vl t + 5z£\t) 
x R x (vxt + 5z\ 1] , t)R 2 (d + v 2 t + 5zf ] , t)} 
= (r\ (vxt + 5z[ 1} , t)Rx (vxt + 5zf- ] , t) \ 

x e i<&(d-&vt+&zf ) -&z^ ) ,t)-i<b(d-Avt+8zf ) -5z\ 1 \t) 



= /r\ (vxt + Sz[ 1} , t)R x (vxt + Szl 1 ] , t) \ 



x e 



i(01-02)(l-j)(fe-i) 



(C13) 



Appendix D: A specific decoherence model 

The following model is inspired by the decoherence 
model for time-dependent Kerr nonlinearities of Ref. [2Tj 
and restricted to the case Vx(z) = V 2 (z) = V(z). The 
reservoir that produces the decoherence is composed out 
of harmonic oscillators in their ground state. The deco- 
herence operators are given by 



rn i (z,t) = M(z,t) + M^{z,t), 

-\M(z,t) + \M^(z,t), 



m 2 (z,t) 
M(z,t) 



(Dl) 
(D2) 



r/Av 
4ttviv 2 



duj (iB(z, u)e~^ 1 + C{z, uj)e luJt ), 

(D3) 



with 



[B(z,u),bHz',u')] = [C(z,w),C\z',u>')} 

= V(z- z')5(u-d). (D4) 



The operators rrii fulfill the commutation relations ( 26 1 



Using the Baker-Campbell-Hausdorff equation to sepa- 
rate the annihilation part M(z, t) and the creation part 
Mt( 2 ,t) in the exponent Eq. (37 1 can be proven. With 



the same method a straightforward but tedious calcula- 
tion yields 



c t = (R t {z,t)) 2 , 
AW = ^oi.io = e5(^-«cic 2 . 



(D5) 
(D6) 
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